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To recover the long-time behavior and the statistics of molecular trajectories from the large 
number (tens of thousands) of their short fragments, obtained by super-resolution methods 
at the single molecule level inEi, data analysis based on a stochastic model of their non¬ 
equilibrium motion is required. Recently, we characterized the local biophysical properties 
underlying receptor motion based on coarse-grained long-range interactions, corresponding 
to attracting potential wells of large sizes |3], as predicted theoretically in d [5] (see also [6] 
for corrals). The purpose of this letter is to point out what was done and not done in d 
this subject. 

Parameter estimation 

The new paper of Masson et ah [7] presents a comparable study of glycine receptors at 
inhibitory synapses. It is our aim here to clarify some issues that arise from the analysis 
and simulations presented there. In particular, [7] asserts that “the method used in Hoze et 
al. (9), also based on a combination of high-density single-molecule imaging and statistical 
inference, evaluates the diffusion and drift by computing the maximal likelihood estimation 
in a mesh sguare as described in Turkcan et al. (7).” This statement may be misleading, 
because the analysis of [3] does not concern “maximal likelihood estimation.” The analysis of 
[3] is based on the assumption that the molecular motion is a diffusion process, whose local 
drift and diffusion coefficients are estimated from local conditional moments of the available 
many fragments of the trajectories. The estimates in |3] are the standard non-parametric 
empirical estimates of the drift and diffusion coefficients known for decades ii EiEni El]. 
Thus the analysis presented in Tiirkcan et al jT2], which was actually published after [3] has 
been sent for publication, does not directly concern the issue at hand. 

The physical model of a receptor motion on a homogenous surface is the overdamped limit 
of the Saffman-Delbriick-Langevin model [T3|, [13]. Specihcally, the diffusion of a receptor 
embedded in a membrane surface is generated by a diffusion coefficient D and a held of force 
F{X), 

X = hW+V^ti,, (1) 

7 

where is a vector of independent standard d-correlated Gaussian white noises and 7 is 
the dynamical viscosity m- However, the trajectory fragments are not collected on the 
microscopic time-scale of ([T]), but rather on a coarser time-scale of the recording apparatus 
that coarse-grains short events that are due to crowding organization of a variety of obstacles. 
Thus (EJ is coarse-grained on a coarser spatiotemporal scale into an effective stochastic 
equation as (see [16]) 

X = b{X) + V2Be{X)w, (2) 
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with the empirical drift held b{X) and diffusion matrix Bg{X), where the effective diffusion 
tensor De(X) = ^Be{X)B'^ (X) can be expressed in terms of the more microscopic diffusion 
coefficient D and in terms of the density and geometry of obstacles. The observed effective 
diffusion tensor needs not be isotropic and can be state dependent, whereas the friction 
coefficient 7 remains constant. Obviously, the impenetrable obstacles that slow the effective 
diffusion down affect neither the microscopic physical properties of the diffusing particle nor 
those of the membrane. Note that the effective held b{X) may have no potential (see |10j). 

Theoretically, the coefficients of ([T]) can be statistically estimated at each point X = 
of the membrane from an inhnitely large sample of its trajectories in the neigh¬ 
borhood of the point X at time t. Specihcally, setting AX{t) = X(t + At) — X(t), the 
inhnitesimal coefficients are given by the limits of the expectations | 10 ] . m 


b{X) 


E[AX(f) |X(f) = X] 

hm —- - 

At^O At 


DI 


E[AX(f)^AX(t) |X(t) = X] 
At^-o 2At 


(3) 


In practice, the expectations are estimated by hnite sample averages and At is the time- 
resolution of the recording of the trajectories, as described in |3]. Thus the estimates [7] 
seem overly optimistic as far as the relation between the theoretical and the sampled tra¬ 
jectories is concerned. By dehnition, the moment estimation in ([3]) requires small fragments 
of trajectories passing through each point of the membrane surface, which is precisely the 
massive data generated by the sptPALM method BE] on biological samples. The empirical 
estimator dehned by ([3]) has been used over the past 60 years in signal processing and 
recently applied in several cell biology contexts HT]. The empirical estimator used in [16] is 
optimal as its variance is minimal and equals to the Cramer-Rao lower bound, as shown in 
section HIE of mi. using numerical simulations. This approach is clearly different from the 
Bayesian method (see formula 6 in [7] and previous publications). 

A positive point of the manuscript [7] is the attempt to include in the analysis the influ¬ 
ence of the Gaussian instrumental noise on the localization of moving receptors and thus on 
their trajectories. However, the variance term in formula 6 of [7] is inapplicable to sort out 
the biophysical parameters when there is a force or a drift component that can vary signif¬ 
icantly in space. This is indeed the case at potential wells. The additional term, needed in 
formula 6 of [7] appears in the direct estimation of the variance E ((X^+i — X„)^|X„ = X) 
and is proportional to the Laplacian of a potential well AV, multiplied by the variance of 
the localization precision. Finally, the motion blur due to the open camera shutter during 
acquisition can also be included in the analysis, as recently discussed in mi. 

The potential energy 

The fluctuation-dissipation principle invoked in |7] expresses the relation between the lo¬ 
cal friction acting on the diffusing receptor and the diffusion coefficient. This relation is 
microscopic and is a part of the collision model of diffusion. The observed coarse-grained 
dynamics, which depends on obstacle density and may vary spatially, while friction remains 
constant, should not be expected to obey the fluctuation-dissipation principle [18]. Formula 
1 of [7] postulates that the friction coefficient is an average with respect to the obstacle 
density, which leads to a coarse-graining approximation, incompatible with the local inter¬ 
pretation that the potential creates by local molecular interactions (which can extend only to 
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few nanometers). This interpretation is thus inconsistent with the hundreds-of-nanometers 
size of the effective potential well described in Fig. 3a of [7]. As described above, the effec¬ 
tive held b{X) that accounts for local traps [5] may have no potential at all. The relation 
between the local molecular energy and the size and depth of the potential well remains 
unclear. 

Energy maps 

The purpose of estimating the drift term is to clarify whether forces other than diffusion 
contribute to the dynamics. In this respect, although a clear large locally conhning potential 
well was presented in [7], it is not clear what is the meaning of the energy map presented 
in Fig. 2c and f. Potential interactions are local and the force is due to the gradient of the 
potential, thus a large region with given constant energy indicates pure Brownian motion in 
this region, without interactions with any held. Is the energy actually calibrated or extracted 
from data? 

Simulations 

Another unclear point of the paper [7] is the role, goal, and the construction of the sim¬ 
ulations. The statement ’’Fokker-Planck equations can always be approximated by master 
equations” seems to indicate that the observed drift vector and diffusion tensor, which may 
be space and time-dependent, are known in the entire domain for all relevant times, which 
may be quite long. The Fokker-Planck equation (FPE) can be extracted from its global 
empirical solution, which is not contained in the empirical data. Thus the only way to 
extract a FPE from the data is to estimate the coefficients from the data, as described in 
[3]. Simulations of long trajectories of the assumed diffusion process can be run once the 
coefficients have been estimated in the entire domain. The approximations of the solution of 
the FPE proposed in [7] is a path integral that requires entire trajectories, not only disjoint 
fragments. The statement ’’these simulated trajectories [...] characteristics match those of 
the experimental ones” should be clarihed, because theoretically, all diffusion processes may 
have the same trajectories, but with different probabilities. Thus simulated trajectories do 
not ’’match” empirical trajectories. What can be matched are various statistics of the tra¬ 
jectories, though it is not enough to match moments, because statistics, such as hrst passage 
times, which may be rare events, cannot be extracted from short simulations and certainly 
not from short fragments of trajectories. The computation of the MFPT of a given diffusion 
process to a boundary requires either solving the Pontryagin-Andronov-Vitt boundary value 
problem HD] or running Brownian simulations, once the drift and diffusion tensor have been 
estimated. But it cannot be computed from the solution of the Fokker-Planck equation in 
free space, because absorbing boundary conditions are necessary. A new approach based on 
stochastic simulations in the entire dendrite is developed in [T5] . 

Finally, estimating biophysical quantities, such as the hrst and second moments from the 
path-integral approach used in [7] may introduce another computational bias: using entire 
trajectories instead of two consecutive points from the same trajectory to compute these 
moments results in an estimation procedure that smoothes out the values of the diffusion 
coefficient. The diffusion coefficients and energy extracted in [7] are signihcantly lower and 
smoother than the ones presented in [3]. It is not clear whether this is a consequence of the 
method of [7] or it captures the biophysical reality of inhibitory receptor trafficking. 
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